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Abstract. Using the standard Green-Kubo formalism, we determine the shear viscosity T] of a hot and dense Yukawa-Fermi 
gas. In particular, we study the effect of particle and plasmino excitations on thermal properties of the fermionic part of the 
shear viscosity, and explore the effects of thermal corrections to particle masses on bosonic and fermionic shear viscosities, Tj^, 
and T)y. It turns out that the effects of plasminos on Tj^ become negligible with increasing (decreasing) temperature (chemical 
potential). 
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INTRODUCTION 

Exploring the nature of the medium created after relativistic heavy-ion collisions has attracted much attention in recent 
years [1]. It has been shown, that after thermalization [2] this medium undertakes a collective evolution. Theoretically, 
this indicates the dominance of the long wave-length limit of the considered theory. The corresponding thermal 
medium can therefore be described by relativistic hydrodynamics or, in general, relativistic viscous hydrodynamics. 
The latter is characterized by a number of transport coefficients. In the present work, we will focus on the shear 
viscosity, which measures the resistance of a fluid against a lateral flow [3]. In general, there are two different 
approaches for determining transport coefficients: The Kinetic Theory (KT) and the Green-Kubo formalism. The 
former is based on solving the Boltzmann equation in the linear regime in which all disturbances are small enough 
compared to the average values [4], and the latter is based on the theory of linear response. In the framework of KT, 
certain divergences appear in the perturbative computation of transport coefficients. In the massless theories and at 
high temperature, these divergences arise from small-angle scatterings [5], and can be remedied by including Hard 
Thermal Loop (HTL) corrections to the internal propagators of the fermionic modes [6]. In the Green-Kubo approach, 
in contrast, these divergences appear as a resummation of an infinite number of diagrams, which has one and the 
same order of coupling constant [7]. In this way, the conventional perturbative treatment breaks down in the long 
wave-length (hydrodynamics) limit of the considered theory [8]. 

In the present work, we have used the Green-Kubo formalism to determine the dependence of the shear viscosity 
of a hot and dense Yukawa-Fermi gas on temperature (T), chemical potential (/r) as well as on bosonic and fermionic 
masses, mi, and niff Since we are interested on the regimes of low and moderate temperature, where particle masses 
cannot be ignored, we shall not be worried about the breakdown of our perturbative treatment. As aforementioned, the 
latter is caused by the appearance of certain divergences in the chiral limit. Considering the bosonic and fermionic 
masses enables us to include HTL corrections via an appropriate modification of particle masses in the thermal 
medium. The novelty in the present work is considering the contribution of particle and plasmino modes [10] in the 
thermal medium, and studying, in particular, their contributions to the fermionic part of the shear viscosity. We will 
show that these contributions significantly affect the shear viscosity, in particular at low and moderate temperatures 

[9]. 
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FIGURE 1. Dashed and solid lines denote dressed bosonic and fermionic two-point Green functions, respectively. According to 
the arguments presented in the main text, we limit ourselves only to the first two diagrams [9]. 


THE METHOD 


Our approach is similar to what is performed in [11]. To start, let us introduce the Green-Kubo type formula for the 
shear viscosity 


ri = ^Jd^x J 


( 1 ) 


Here, jS* = yjS is the inverse proper temperature, with y = 
momentum tensor, 7’^'', 



j, and is the traceless part of energy- 




( 2 ) 


To determine rj from (1) in a diagrammatic manner, we adopt an appropriate skeleton expansion, which is demon¬ 
strated in Fig. 1. Here, the lines correspond to dressed two-point Green functions (TPGF) of bosons (dashed lines) and 
fermions (solid lines). In (2), Tpa is the energy-momentum tensor of the underlying quantum field theory, and A^'' is 
an operator which project any quantity in hypersurface orthogonal to the direction of fluid velocity . It is defined by 
AP'' = gP'' — uPid', where the metric = diag(+,—, — 

As it turns out [11, 9], Tj is sensitive to low frequency limit of the considered theory. It is given by 




po^o 


(3) 


where np(/7o) is the Fourier transform of the retarded correlator of two traceless operators k^'' from (2). According to 
the above formulation, the shear viscosity tj can be determined by adopting a specific theory, deriving the traceless part 
of its energy-momentum tensor from (2) and inserting it into (3). In the present work, we will select the Yukawa theory 
and will determine the corresponding bosonic and fermionic shear vicosities. The Yukawa theory has the well-known 
Lagrangian 


jSf = \if{i ^ 




The corresponding energy-momentum tensor of this Lagrangian reads 


(4) 


Tpv — T ~ -S ?gpv 


(5) 


Any retarded TPGF of this theory can be related to the corresponding spectral function, p, via the Kallen-Lehmann 
representation 


Gr(p,®„) = J 


dco p{p,co) 
2n CO + i(0„' 


( 6 ) 


Using (6), the spectral function can be shown to be proportional to the imaginary part of the corresponding retarded 
TPGF. Similarly, the real part of each retarded TPGF can be derived using Kramers-Kronig relation. In this way, the 
shear viscosity is expressed in terms of the real and imaginary parts of retarded bosonic and fermionic TPGFs. 

At this stage some remarks are in order: For the bosonic particles, because of certain symmetry properties [12], we 
have in general no ambiguity to determine the spectral function, pp. It is given by [11,9] 


Pb{p,(o) 


_ 4-(orh{p,(Oh) _ 

- eKp, cOb) - Flip, (Ob)]^ + 4(02rl{p, (Ob) ’ 


(7) 






where a)| = E^{p) = ri,{p) = — 2 ^Im[Z^(p)]. Here, is the retarded bosonic TPGF. 

As concerns the fermionic spectral function Pf, however, certain ambiguity occurs. As it turns out, Pf includes both 
particles and antiparticles. In the chiral limit, it is given by [13] 


Pf 


(®-£,)2+r2/2 


h+{k)- 


{co + Ek)^ + rl/2 


h-ik), 


( 8 ) 


where the projector h+ (k) and /z_ (k) correspond to particle and antiparticle contribution to pf, respectively. Moreover, 
is the decay width of both fermions and antifermions. Let us notice that at high temperature, where nif <ffT, there is 
no difference between the decay widths corresponding to fermions and antifermions. However, at moderate T, where 
particle masses turn out to have nontrivial effects on the quantities appearing in Pf in ( 8 ), we have to distinguish 
between the decay widths of fermions and antifermions. One of these effects is the appearance of new degrees of 
freedom, the so called plasminos [14]. Plasminos are similar to antiparticles, in the sense that both plasminos and 
antiparticles have opposite chirality and helicity. Their difference arises in their corresponding dispersion relations, 
once the contributions of T-dependent part of the radiative corrections are taken into account [14]. The expression ( 8 ) 
for Pf turns out to be therefore invalid for nif k,T . Here, the difference between particles and plasminos implies a 
more accurate form of Pf [9], 


P/(P,®) 


_ 2r+(p,a)/) __ 2r-(p,CQ/) _ 

(®-£+(p,a )/))2 + r^(p,®/)^+ (® + £-(p ,®/))2 + r 2 (p,®/)^^’ 


(9) 


where (dj = and 

KOf 

In (9), E±{p,(Of) and r±(p,®^) denote the real and imaginary parts of the fermionic TPGF, respectively. They are 
defined by 


1 f 

£±(p,®/) = ®/±-fr(|±(p,®/)Re[E^(p,®/)]), 

r±(p,®/) = ±^fr(|±(p,®/)Im[E^(p,®/)]). (10) 

Let us notice that in (9), comparing to ( 8 ), particle masses as well as the difference between the decay widths 
corresponding to particles and plasminos are taken into account [9]. Here, the imaginary part of bosonic and fermionic 
TPGT is computed using the real-time formalism (RTF) of finite temperature field theory (TFTF) [16, 17]. In what 
follows, we will present the final results for r;,,r± as well as bosonic and fermionic shear viscosities 77 * and r\f. 


ANALYTICAL RESULTS 

According to the arguments presented in the previous sections, 77 is sensitive to the low frequencies. By decomposing 
(1) into bosonic and fermionic parts, we can calculate 77 * and r\f. The retarded bosonic and fermionic TPGF extracted 
from the first two terms of Fig. 1 read [9] 

nfo(po) = -^^ri'^^^'^PaPpPpPa j ^^^^^Pfo(p,®i)p*(p,®2K(®l)n*(®2)We(®i2,po), (H) 


Here, ri‘^PP^ = — ^AP'’AP^, 0)12 = COi + (O 2 , CO 12 = COi — (O 2 , and 


We((0,po) 


1 1 
Po + is — CO po + is + CO 
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FIGURE 2. (color online), (a) The T dependence of U;, for co/, = 300 MeV and (b) the /l dependence of F/, for co/, = 300 
MeV. The dashed lines correspond to T;, including the constant mass contributions of bosons, while the solid lines correspond to 
the same quantity including thermal corrections of bosonic masses described by (17). The red, gray and blue lines correspond to 
ra® = 100,150,300 MeV and ra® = 5 MeV. Here, the Yukawa coupling g = 0.5 is used [9]. 


Moreover, ni,^f{co) are bosonic {b) and fermionic (/) distribution functions, given by 

rihico) = -, nf=—;i -. ( 12 ) 

The bosonic and fermionic spectral functions, Pi(p,a)) and Pf{p,co), are defined in (7) and (9), respectively. By 
inserting these relations into (3), and after some straightforward computations, the bosonic and fermionic shear 
viscosity rjh and rjf read (see [9] for more details) 

Tib = 

nr = 

In the fermionic channel, i = ± refer to the particle (s = +) and plasmino (s = —) excitations, respectively. In addition, 
r* are defined by = r_)_ ± r_. We are working in a frame, in which quantum fluctuations do not imply any sensible 
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+ ^(r±). 


(13) 




FIGURE 3. (color online). The T dependence of rjt, is plotted for p. = 120 MeV for (a) (T, p)-independent <^o = 40,80 and (b) 
(T, p)-dependent = 40,80 [9]. 






























a) cof = 300 MeV, // = 150 MeV 


b) ojf = 300 MeV, = 150 MeV 
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FIGURE 4. (color online). The T dependence of (a) r+, (b) r_ is plotted for constant cOf = 300 MeV and /i = 150 MeV. Dashed 
lines correspond to r±, including (T, ^(-independent = 60,90,120, while the solid lines correspond to the same quantities 
including the HTL corrections to the bosonic and fermionic masses [9]. 



effects on dispersion relations. We therefore adopt the approximations Ei, Ri O/, and Ef r; (Of, and compute F;, as well 
as r± using the RTF of FTFT up to hrst order correction of a perturbative expansion in the orders of the Yukawa 
coupling constant g. Analytical results for the bosonic and fermionic decay widths are given by 


r^Cp,®) 




-cosh(T/) +cosh^ + 
■cosh(T/) +cosh^ M — 


4)(l-7|)y 


(14) 


and 


r+(p,®) = 


g27’7/(^^-4) 


32n 
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FIGURE 5. (color online). The p dependence of (a) r+, (b) F- is plotted for constant (Of = 300 MeV and T = 150 MeV. Dashed 
lines correspond to r±, including (T, p (-independent = 60,90,120, while the solid lines correspond to the same quantities 
including the HTL corrections to the bosonic and fermionic masses [9]. 
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FIGURE 6. (color online). The T dependence of rjf is plotted for fl = 120 MeV and for (a) (T,)i)-independent <^o = 40,80 as 
well as (b) (T,jU)-dependent = 40,80 [9]. 


Moreover, we have r_ = r_|_ — with 


r? = - 
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1 -H- 2 — - 2 — 


• (16) 


In the above relations, Jb = 'xib/fOb^Yf = ^ = mhimf as well as Xf = = a)/,j3 and Kf = COfp. We 

have also used the definitions 


s± = ± _ 4)(1 _ ^ 2 )] / 2 , and Y± = jf/( 7 /± l)/2. 

To consider HTL corrections for bosonic and fermionic masses, we have introduced = mi,{T, jj.)/mf{T, jj.) with 
mh(T,n) = m^ + nY^'{T,n) and mf{T,jj.) = + m'^(Y,/r). Here, the superscript 0 corresponds to the T- and jx- 
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FIGURE 7. (color online). The T dependence of A is plotted for p = 120 MeV and ^o = = 40. 













































independent bosonic and fermionic masses, whose thermal corrections are given by 



Let us notice at this stage that the coupling g does not receive any HTL corrections in the Yukawa theory. In this 
section, we have presented the basic analytical expressions, which are necessary for calculating the decay widths 
and r± as well as the shear viscosities ribjilf terms of free parameters of the theory. The T and jj. dependence of 
the above quantities will be demonstrated in the next section. 


NUMERICAL RESULTS 

The T and /r dependence of is demonstrated in Fig. 2. As it turns out, by increasing T and jj., F;, decreases. 
This reflects the fact that the bosonic particles have a larger mean free path in dense or hot medium. By considering 
thermal corrections to the bosonic mass, the mean free path of bosons increases. The T and jU dependence of rjb can 
be determined by inserting the expression for F^ into the (13). In Fig. 3, the T dependence of rjb is demonstrated for 
jj. = 120 MeV and for T- and /r-independent and dependent bosonic to fermionic mass ratios, and . The fact that 
rjb increases with T demonstrates the fact that the mean free path (time) of bosons increases in a hot medium. This is 
consistent with the arguments based on relaxation time approximation. By including HTL corrections to the particle 
masses at fixed /r, rjb decreases for each fixed T and jj.. 

In Figs. 4 and 5, the T and /r dependence of the decay widths corresponding to particle and plasmino modes are 
plotted. The T and /r dependence of r+ is somehow similar to that corresponding to F^. In contrast, the T and /r 
dependence of r_ seems to be different, in the sense that by increasing T and jj., r_ increases. This indicates that at 
high temperature or in dense medium the plasmino modes wash out earlier than the particle modes. In particular, this 
means that by increasing T and jj., plasmino modes decay into the normal particle modes, and hence the normal particle 

f 

lifetime increases with increasing temperature or chemical potential. By plugging F^ into (13), the T dependence of 
rjf can be determined for different = 40,80 and = 40,80 and fixed jj.. The results of T]^ are demonstrated in 
Fig. 6. Unlike the bosonic case, the fermionic shear viscosity is a decreasing function in terms of T. This denotes that 
the mean free path (time) of fermionic particles decrease with increasing temperature. Moreover, by including HTL 
corrections for particle masses, rjf increases for each fixed T and jj., in contrast to the bosonic case. 

Our main result is presented in Fig. 7, where the T dependence of a certain quantity A is presented. Here, A is 
defined by the difference between rjf in terms of r+ = r_ and r+ ^ r_, i.e., 

A = 77y[r+ = r_]-77y[r+^r_]. (i8) 

Figure 7 illustrates the fact that by increasing temperature, the effect of plasmino modes on rjf become negligible, and 
the approximation r+ « r_ is therefore reliable only in the limit of high temperature [13, 18]. 


CONCLUDING REMARKS 

The shear viscosity is a transport coefficient which quantifies the long wave-length response of a medium to momentum 
anisotropies. In the present work, we calculated the bosonic and fermionic parts of the shear viscosity of a hot and dense 
Yukawa-Fermi gas in the leading ^(F^*) order. Since our theory is massive, the small-angle scattering divergences 
are irrelevant, and a perturbative computation of transport coefficients is reliable. We showed that by increasing the 
temperature, Vjb increases. In contrast, Vjf decreases with increasing temperature. We also took the HTL corrections 
to particle masses into account and investigated their effects on Vjb as well as Vjf. We eventually showed that the 
assumption F+ « F_ is only justified at high temperature. Here, F+ and F_ are the decay widths corresponding to 
particle and plasmino excitations. 
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